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In the path integral approach to false vacuum decay with the effect of gravity, there is an un¬ 
solved problem, called the negative mode problem. We show that the appearance of a supercritical 
supercurvature mode in the one-bubble open inflation scenario is equivalent to the existence of a 
negative mode around the Euclidean bounce solution. Supercritical supercurvature modes are those 
whose mode functions diverge exponentially for large spatial radius on the time constant hyper¬ 
surface of the open universe. Then we propose a conjecture that there should be “no supercritical 
supercurvature mode”. For a class of models that contains a wide variety of tunneling potentials, 
this conjecture is shown to be correct. 
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I. INTRODUCTION 


The Euclidean path integral approach has been used to investigate the true vacuum bubble nucleation through 
quantum tunneling [QH]. In the lowest WKB approximation, the quantum tunneling is described by a bounce solution. 
The bounce solution is a solution of the Euclidean equation of motion, which connects the configurations before and 
after tunneling. The bounce solution that takes account of the gravitational effect was found by Coleman and De 
Luccia ||]. 

The decay rate per unit volume and per unit time interval, T, is given by the formula ^ 


r = 


\K\e 


-T 


(bounce) f^(trivial) s. 


( 1 . 1 ) 


where is classical Euclidean action for the bounce solution and fgj. trivial solution 

that stays at the false vacuum. In the path integral approach, the prefactor \K\ is evaluated by the gaussian integral 
over fluctuations around the background bounce solution. In a standard system which does not take account of 
gravity, there is one perturbation mode in which direction the action decreases. It is called a negative mode. For this 
mode, the gaussian integral is not well-dehned. To make the integral finite, the integration path should be deformed 
on the complex plane. Consequently, one imaginary factor, i, appears in K. In the Euclidean path integral approach 
to tunneling, this imaginary unit plays a crucial role to interpret T as the decay rate. 

However, in the case when gravity is taken into account, the situation changes drastically. Since there are gauge 
degrees of freedom, we have various possibilities in choosing variables to describe the physical degrees of freedom. If 
we choose inappropriate variables, the equation that determines the fluctuation mode can become singular. In Ref. 
0 , we have shown that it is impossible to obtain a well-behaved reduced action as long as we stick to variables that 
appear in the original Lagrangian in the second order formalism. In order to hnd variables that can lead us to a 
well-behaved reduced action, it was necessary to resort to the Hamiltonian formalism. In the Hamiltonian formalism, 
conjugate momenta are introduced, and hence wider varieties of choice of variables are allowed. In Ref. Q, we derived 
a well-behaved reduced action for fluctuations around the 0(4)-symmetric bounce solution. We found the action 
for fluctuations which conserve the 0(4)-symmetry has an unusual signature. Namely the kinetic term is negative 
definite. To deal with this action, we proposed a prescription analogous to the conformal rotation [^. Then, from the 
path integral measure, there arises one imaginary unit i. This suggests there should be no negative mode in the final 
form of the reduced action, since the factor i in K has been already taken care of by the above mentioned prescription. 
Therefore, we proposed the “no negative mode conjecture” [Q. 
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On the other hand, in recent years the process of false vacuum decay with the effect of gravity has been studied 
extensively in the context of the one-bubble open inflation scenario. In this scenario, an open universe is created inside 
a nucleated bubble. In one-bubble open inflation, one of the most important issues is to calculate the spectrum of 
quantum fluctuations after the bubble nucleation because it determines the spectrum of cosmological perturbations. 
By comparing the predicted spectrum with the observed one, we can test a model of one-bubble open inflation. 

Fluctuations in an open universe can be decomposed by using spatial harmonics on the unit 3-dimensional hyperbolic 
space. We denote the eigenvalue of a spatial harmonic by —{p^ + l). The spatial harmonics with positive are square- 
integrable functions on a time constant hypersurface in an open universe in the sense that they can be normalized 
by using the Dirac delta function. As a result the spectrum is continuous for > 0. On the other hand, the spatial 
harmonics are no longer square-integrable for p^ < 0. However, since a time constant hypersurface in an open universe 
is not a Cauchy surface, this divergence does not directly exclude such modes. By considering the normalization of 
perturbation modes on a Cauchy surface, we find that the spectrum for p^ < 0 becomes discrete. These modes are 
called supercurvature modes since they give rise to correlation on scales greater than the spatial curvature scale . 

There are two classes of supercurvature modes, which we call supercritical and subcritical modes. The precise 
definition will be presented later. Here we mention that supercritical modes have smaller values (i.e., larger absolute 
values) of than subcritical modes, and their mode functions diverge exponentially for large spatial radius in the open 
universe. We shall show that the existence of a supercritical supercurvature mode is equivalent to the existence of a 
negative mode. Thus the “no negative mode conjecture” presented in Ref. Q| can be restated as the “no supercritical 
supercurvature mode conjecture”. Note that there exists no supercurvature mode for a model with sufhciently thin 
bubble wall |^, hence there is no supercritical mode in such a model. The purpose of this paper is to examine this 
conjecture for a wider class of models. In particular, a class of models we consider naturally allows thick bubble wall 
solutions. Recently the negative mode problem has been discussed by Lavrelashvili ||] in the context of the singular 
Hawking-Turok instanton ]lC| ] . In this paper, however, we exclude the possibility of the Hawking-Turok instanton and 
focus on regular bounce solutions. 

In section 2, we show the equivalence between the existence of a negative mode in the reduced Euclidean action 
and that of a supercritical supercurvature mode in the spectrum of cosmological perturbations. In section 3, we give 
a method to construct potential models which allow an analytical treatment. By using this method, we give a set 
of models which seemingly violate our conjecture. In section 4, for such models that have a bounce solution with a 
negative mode, we show that there should be another bounce solution that has a smaller value of the action and has 
no negative mode. Conclusion is given in section 5. 


II. NEGATIVE MODE AND SUPERCRITICAL SUPERCURVATURE MODE 

We consider a system composed of a real scalar field, $, coupled with the Einstein gravity. The Euclidean action 
is given by 


Se= / 


2 k 


I 


+ U($) 


( 2 . 1 ) 


The potential of the scalar field is assumed to have the form as shown in Fig. IH and initially the field is assumed 
to be trapped in the false vacuum. As mentioned in Introduction, the bounce smution is a Euclidean solution that 
connects the configurations before and after tunneling. In the present case, the geometry before false vacuum decay 
is given by a de Sitter space. After tunneling, there appears a true vacuum bubble in the false vacuum sea. This 
bounce solution is obtained by Coleman and De Luccia 0 under the assumption of the 0(4)-symmetry: 


= a^ip) -I- dx^ + sin^ x(iD(2)} > = ^iv)- 

The Euclidean equations of motion are 
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(2.3) 

(2.4) 

(2.5) 


where a prime represents the differentiation with respect to p and TL := a'/a. Requiring the regularity of the bounce 
solution, the boundary condition is determined as 


(j)'(X e a oc e (ry ^ Too). 
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( 2 . 6 ) 
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By choosing the initial value of at 77 ^ —00 appropriately, we obtain a solution which satisfies the required boundary 
condition. For definiteness, we choose </> to be in the false vacuum side for 77 —> — 00 . 



FIG. 1. A typical shape of a scalar field potential under consideration. 


As noted before, the number of I’s in the prefactor K has a crucial meaning in the path integral approach. To 
evaluate this number, we need to derive the reduced action for fluctuations around the bounce solution. The fluctua¬ 
tions can be expanded in terms of the spherical harmonics on the unit 3-sphere, say Ynim with n = 0, 1, 2, • • •, which 
satisfy -I- n{n + 2)]Ynem = 0. After appropriate gauge fixing, we obtain p|,[ll|] 


(n + 3)(n-l) r 


n£m 


^nim 


dq 


.nim 


drj 


+ iU + {n + ?,){n-l))\q 


nim 1 2 


(2.7) 


where is the coefficient of the harmonic expansion of a gauge-invariant variable g, and its conjugate. In 

our original derivation [Q, we used the variable h := q/a and its conjugate momentum. The variable g is equal to the 
Euclidean version of q introduced in Ref. Here, the potential U is given by 


U = ^(^'2 J^^r2 W=^. 

2 (j) 


( 2 . 8 ) 


Physically, the variable g represents the curvature perturbation in the Newton gauge when analytically continued to 
the open universe inside the bubble. The explicit form of the perturbation in this gauge is 


ds‘‘ = 


1-H 


K(j)'q 


a 


dr]^+[l- 


K(j)'q 


-f sin^ Xd^l2)) 


where is the perturbation of the scalar held. Recall that the prefactor K is evaluated as |Q 

K = J [dTrdgJe”^' 


(2.9) 


( 2 . 10 ) 


Looking at Eq. (^), we hnd that the action disappears for n = 1. This just reflects the fact that the n = 1 modes 
are pure gauge and there is no physical degree of freedom in them. This is very different from the case of false vacuum 
decay in the hat spacetime, in which there is a zero mode in the n = 1 modes. This zero mode describes spacetime 
translation modes and its existence implies the existence of a unique negative mode in the n = 0 (i.e., 0(4) symmetric) 
modes ||^. In the present case of false vacuum decay with gravity, the coefficient in front of | 7 r "-^™|2 becomes negative 
for n = 0. Therefore if we try to perform the integration with respect to tt, we hnd that the gaussian integral does not 
converge for n = 0. To resolve this difficulty, we proposed in 0 a prescription analogous to the conformal rotation. 
By changing the variables as tt ^ —in, q —>■ iq, the above path integral becomes well-dehned. To carry out the path 
integral, the variables must be discretized. Then the numbers of p and g integrals will differ by 1. Therefore, this 
change of variables will produce one imaginary unit i for the prefactor K. This is in contrast to the case of false 
vacuum decay in the flat spacetime in which the imaginary unit i in K arises from the negative mode. Note that 
there is no negative mode for n > 2 H. 


Let us examine the n = 0 case in more detail. 
Shrodinger-type equation. 


If we know the spectrum of eigenvalues, \j, of the following 
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d?7" 


+ C/-3 


_ /2 

where qj is the j-th eigenfunction, the contribution to K from n = 0 modes, will be given by ^ (—i) ' , where the 

factor (—i) comes from the rotation of the variables as explained above. If there is a mode with a negative eigenvalue, 
there arises another imaginary factor, and hence the prefactor K becomes real. If so, this bounce solution will not 
contribute to the tunneling process in the path integral approach. On the other hand, there is another approach to 
describe the tunneling. That is the approach to construct a WKB wave function . In this approach, any bounce 
solution that connects false and true vacua will contribute to the tunneling, putting aside the issue of which bounce 
dominates. Thus, in order that the consistency is kept between the path integral approach and the wave function 
approach, there should be no negative mode for the Coleman-De Luccia bounce solution, provided it gives the smallest 
Euclidean action among the non-trivial bounce solutions. This is the conjecture proposed in Ref. 





FIG. 2. Conformal diagram of a de Sitter-like space with a bubble wall. The Cauchy surface exists in region C but not in 
the open universe inside bubble, the region marked as R. 


Now, let us turn to the problem of the perturbation spectrum in the context of the one-bubble open inflation 
scenario. To study it, we need to quantize the perturbation field on the non-trivial background that appears after 
tunneling. For this purpose, the reduced action for the perturbation field must be calculated. In doing so, the time 
coordinate must be chosen so that the time constant hypersurface is a Cauchy surface. A convenient choice is to use 
the coordinate that is obtained by the analytic continuation of Xc = *(x ~ '^/2)- Then, the coordinates (xo ?7, ^( 2 )) 
span the region C of Fig. ||. In region C, the metric takes the form 


ds^ = a?{r]) — dxl + cosh^ Xc ■ 

It should be remembered that rj is not a time coordinate there. The reduced 
1 


H 9 / 


cosh^Xc^O^-g^^O 
OXc OXc 


+ 1 ) + 


action in this region is 

(0-3)cosh2xc}g'™ 


( 2 . 12 ) 

obtained as 

(2.13) 


where O is defined by 


d 



(2.14) 


* Rigorously speaking, there can be even number of negative modes, depending on the choice of the canonical variables. 
However, for the present choice of the variables (tt, q), we conjecture that there is no negative mode. 
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and U is the same as the one given in Eq.(2.8). In this expression, the ry-dependence appears only through the operator 
O. Thus the perturbation can be expanded by using eigenfunctions of O: The equation 

that determines eigenvalues and eigenfunctions is 


dr]^ 


U 


qP = (p2 + 4)qP. 


(2.15) 


Analyzing the asymptotic behavior of (j), we find {7^4 for rj =Foo. Hence the spectrum is continuous for > 0. 
On the other hand, the spectrum for < 0 is discrete, which are called supercurvature modes. Here we call modes 
with p2 < _i supercritical supercurvature, while we call modes with 0 > p^ > — 1 subcritical supercurvature. When 
analytically continued to the op en un iverse by \ r = Xc + the supercritical modes diverge exponentially for 

Xr oo. By comparing Eqs. (2.11) and ( ^.15 ), the existence of a negative mode around the bounce solution is 
manifestly equivalent to that of a supercritical supercurvature mode in the spectrum of quantum fluctuations in the 
open universe. Thus the conjecture that we stated above can be restated as the “no supercritical supercurvature 
mode conjecture”. 


III. MODEL CONSTRUCTION 


At hrst sight, there seems no reason to deny the possibility of a potential model that allows a bounce solution with 
a supercritical supercurvature mode. In fact, there is a method to construct such models, as described in a separate 
paper 0. However, as mentioned there, this method does not guarantee that thus obtained bounce solution gives 
the smallest Euclidean action among non-trivial solutions. Hence, given a bounce solution that allows a supercritical 
supercurvature mode, if one can prove that its action is not the smallest, our conjecture remains intact. 

In this section, we review the method to construct a potential model and a bounce solution, and give an example in 
which a supercritical supercurvature mode appears when a model parameter is varied. We then construct a potential 
model that allows a continuous series of bounce solutions under the weak gravitational backreaction approximation, 
which will be used for the investigation of the negative mode problem in the next section. 

Our method for model construction is as follows. We first give (j)'■ It can be specified rather arbitrarily except for 
the conditions that it vanishes as for r ? Too and it is positive (or negative) dehnite so that (j) is a, monotonic 

function of rj. Then we solve the equations ( p^ ) to find H and a. For rj Too, H is solved as 


n-i 
n + i 


C^e 


2ri 


Thus the required boundary condition, H ±1 for ij —> Too, is automatically satished. Once we know 77, we can 
determine a by integrating 77. Here an integration constant appears, w hich determines the absolute value of the 
potential, say the potential at the false vacuum Vp- Finally, using Eq. (^), the potential V is determined as a 
function of rj as 


vWv)) 


6 — 677^ T 

2Ka2 


(3.1) 


By assumption, since (/) is a monotonic function of rj, by inversely expressing 77 in terms of cj), the shape of the potential 
is uniquely determined. 

Next we turn to the eigenvalue equation (2.11). If there is a negative mode, the solution with A = 0 must have, at 
least, one node. Hence, it is sufficient to examine the equation for A = 0, which is written as 


dr]"^ 


^ ^2 _ 3 


9 = 0. 


(3.2) 


If the second term, K(j)'‘^/2 is large enough, the potential will stay positive and hence q will have no node. Since 
the change in the overall amplitude of (j)' does not alter the function W, we can construct a model in which the 
second term is negligibly small without changing the other terms. Hence, we neglect this positive definite term for 
the moment. 

From the boundary condition, we find that W behaves as 

W ^± 2 , ( 77 ^ too). 


Except for this boundary condition, W can be arbitrarily chosen. Then, by choosing W to stay close to zero for a 
sufficiently long interval of 77 , the potential term, —W' T — 3, will have a wide negative-valued region, and the 
solution for q will have a node. In contrast, if we consider a sufficiently thin wall, the W term will dominate the 
potential when \W\ is small and W < 0. Then, the potential will stay positive, and hence q has no node ||]. 
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For instance, let us consider a set of models in which W is given by 

W = — 2 tanh ((77 — C)/A), 


(3.3) 


where the model parameters ( and A represent the wall location and its thickness, respectively. This is realized by 
setting 




(cosh(( 77 - C)/A)) 


2A ’ 


where /i is a constant. Then, neglecting the K(/)'^/2 term, Eq. (|3.2|) becomes 






1 


cosh ((r? - C)/A) 


+ 1 


g = 0. 


For A = 1, this equation has a regular nodeless solution. 


1 


9 = 


cosh(77 — C) 


(3.4) 


(3.5) 


(3.6) 


Thus the lowest eigenvalue is zero for A = 1. Since the potential term in Eq. ( |3.5| ) is a monotonically decreasing 
function of A, there will be a negative mode for A > 1 . Now let us turn on the K(/)'^/2 term while keeping W the 
same. As we increase fi, the critical value of A for the existence of a negative mode will increase. Eor > f, a 
negative mode will disappear completely for any value of A. 

To make the problem tractable, let us co nside r the weak backreaction limit, in which the geometry can be approxi¬ 
mated by a pure de Sitter space. Erom Eq. (^), we see that this will be the case if Kcj)'^ <C 1. In the above model, we 
have seen that a negative mode will disappear if the k (^'^/2 term gives a large contribution (of order unity or greater) 
to the potential. This will be true for most conceivable potential models. Hence we may restrict our attention to 
a class of models that automatically satisfy the weak backreaction condition when investigating the existence of a 
negative mode. 

In the first non-trivial order of the weak backreaction approximation, we find that the negative mode problem can 
be analyzed on a fixed de Sitter background geometry. Let us explicitly show this fact. 

Eirst we examine the background field equations. Neglecting relative errors of 0{k4>'^), Eq. (2.5) can be readily 
solved as 


1 


at = 


H cosh 77 ’ 


Tib = — tanh 77 , 


(3.7) 


where H is defined by := {n/?>)Vtop where Vtop is the value of V at the top of the potential barrier. It is worthwhile 
to note that the weak backreaction approximation does not imply small potential difference between the false and true 
vacua. The vacuum energy inside the bubble can be much smaller than Vp as long as the bubble radius is sufficiently 
smaller than . However, the weak backreaction approximation does imply Vtop — Vf ^ Vf- Hence we may replace 
Vtop with Vp in the definition of H. H ere we have chosen Vtop for later convenience. Then substituting the lowest 
order solutions of a and Ji to Eq. (^), the scalar field equation becomes 


+ 27fb<))'- ag —((5E) = 0, 
a 0 


(3.8) 


where bV = V — Vt or,. T his gives the lowest order solution of </>. Introducing the variable N := logo, and expanding 


it as Nb + SN, Eq. (|2.5| ) gives the equation for the first order correction to a: 

cosh^ r] dr] \ dr] J 2 ^ 


(3.9) 


Since ~ for 77 ^ ± 00 , the above equation can be integrated and we find SN = 0{K(f>'^) everywhere. If we 

substitute this SN back into Eq. (2.3), we can calculate the correction of 0{K(j)'^) to (f>. 

Next we evaluate the action under the weak backreaction approximation. The Euclidean action for an 0(4)- 
symmetric configuration is given by 


Sp = S'g[a] -I- S'm[Q 


S'g[a] ;= 2 ^^ 


dr] 


^ /2 ^ 2 I 4t/ 

- a - a + a Vt 

K K 
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Sm[a,4’] := 27 ^^ y 


dr] 
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'2 


+ a^5V{(j)) 


(3.10) 

(3.11) 
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Inserting + 6a to the gravitational part of the action, where 5a = abSN = we find 


SgK + H = Sgiat] (1 + 0 ((ac0'")')) , (3.12) 

which follows from the fact that at is a solution to the equation obtained from the action 5*9[a]. Hence, to the first 
order of we have 


S E — Sg [u^] + Sm {ab , (6] , 


(3.13) 


where Sm[ab,4'] is 0{k4>'‘^) relative to <S'g[a{,]. Thus the action evaluated on a fixed de Sitter background is correct 
to the first non-trivial order. Since the background geometry is fixed, we focus on the matter part of the action and 
denote it by in the following. 

Fin ally, we consider 0(4)-symmetric fluctu atio ns around the bounce solution. We rewrite the eigenvalue equation 
( ^.11 ) by using the variable ip introduced in (^). We find 




2H- 




f(/)'2-A-3 




+ A - - a^d^V + 




f (/)'2 - A - 3 


+ = 0 , 


(3.14) 


where d^V = <PV{(f)/dcj)^. Unless A is close to —3, this equation reduces in the weak backreaction limit to 




+ 2Hbp' — ald^{SV)p = —Xp, 


(3.15) 


Here we note that the term a^d'^^SV) can be as large as Ti-l even in the weak backreaction limit. The above equation 
is the same as the equation for 0(4)-symmetric scalar field fluctuations on a fixed de Sitter background, i.e., without 
taking account of the metric perturbation. 

In the above, we have shown that it is sufficient to consider the problem on a fixed de Sitter background in the first 
non-trivial order of the weak backreaction approximation. Now let us return to the model with W = —2 tanh (?7 — C). 
We assume C > 0 since the 3-sphere of the maximum radius should be on the false vacuum side in order for the bounce 
solution to describe false vacuum decay. As mentioned there, A = 1 is the critical case in this r nodel in the weak 
backreaction limit. This was sho wn b y the fact that there is a regular, nodeless solution of Eq. (Enl) with A = 0 
for A = 1, which is given by Eq. (3.6), and that the po tent ial C/ is a monotonically decreasing function of A. The 
solution for p that corresponds to the nodeless solution (|3.6[) for q is given by 


p oc 


sinh (?7 — C,) 
cosh ^(?7 — C) 


(3.16) 


This solution has a node. In fact, Eq. (t5.I5|) has a solution p = (j)'/ab with a lower (hence negative) eigenvalue A = —3. 
If the gravitational effect were completely neglected, this would be the nega tive mode that gives rise to the factor i 
of the path integral. However, as we have noted, the approximate eq uatio n (3.15) is no longer valid for A « —3 even 
in the weak backreaction limit. Since the unapproximated equation (f5.14|) for p is not a Schrodinger-type equation, 
the corresponde nce b etween the number of nodes and the increasing order of eigenvalues is lost. Thus the existence 
of the solution (3.16) with a node at A = 0 does not imply the existence of a negative mode. 

For A = 1, the background solution for (p is given by 


4> = (pc '■= tanh (?7 - C) + C 2 ), 


(3.17) 


where is a constant which determines the typical scale of (p and ci and C 2 are integration constants to be determined 
later. The potential SV is reconstructed as 


5V := 


1 


6V 




dr] d{6V)p/ 


H^cp- 




= cl 


= —c 


dv — ^ —- (— 2 tanh (?7 — C) ~ 2 tanh? 7 ) 


cosh ‘^(?7 — (/) 
2 smh^(27?- C) 

^ cosh^(r? — C) 
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^2 


(sinh ( tanh (77 — C) + cosh — 1 


n 2 


(3.18) 


2 sinh C L 

where we have set the integration constant so that 5V = 0 at the top of the potential. Then by choosing Ci and C 2 as 


1 . cosh C + 1 

Ci = -sinhC, C2 = - - -, 

we find that all the solutions parametrized by C, are solutions for the same potential given by|] 

2H^ 


SV = -- 




(3.19) 


(3.20) 


Note that the eigenfunction for the A = 0 mode is also given by taking the derivative of the background solution with 
respect to C; oc dcjj^/d^ oc sinh (77 — C)/ cosh^fr; — C)- For C = 0) the scalar field stays still on top of the potential 
barrier. This is a Hawking-Moss type solution |^. It is easy to see that the action S'm[<('c=o] vanishes. Since cj)(^ is a 
continuous series of solutions, the action must be the same for all the values of C- This implies S'm[</'c] ~ ^ for any 
This can be also checked by direct substitution. 


IV. PERTURBATION ANALYSIS 


Let us denote the matter part of the action, corresponding to the critical case of the model discussed in the 

previous section by 

S'(o)[^]:=y’ dr] . (4.1) 

We consider a set of models whose action is given by 

[<),], (4.2) 

where iSF) [(jj] is a small perturbation of order e ^ 1. We need not specify the explicit form of [(/)]. In the following, 
we neglect the O(e^) terms. 

For convenience, we introduce the following notation for variations of the action: 


dr] 


(50 ■ (50(7;) 

>;2 


SS 


Hid) 




^[ 0 ; • 01 , 02 ] ■■= J dr]idT]2 

( 50 (? 7 ) 50 ( 7 ;') 




SHm)mri2) 

0 (^ 0 , 


0l(?7l)02(?72), 


(4.3) 


and analogous expressions for higher variations. 

Let us first write down useful formulas for our discussion below. Since 0,^ is a solution to the unperturbed equation 
of motion for any value of 0 , successive differentiations of the equation of motion with respect to 0 give 


(5S'(o) 

(50 

( 52 ^( 0 ) 

502 

535 .( 0 ) 

503 


0 c; * 


= 0 , 


' 00 ^ 


0 c;*. 


= 0, 

00^ d0(J 


525 ( 0 ) 

502 


0c;* 


02 0 c 


002 


= 0 , 


(4.4) 

(4.5) 

(4.6) 


^Jaume Garriga pointed out to us that this series of instantons can be obtained by a conformal transformation of the Fubini 
instanton in flat space |l4]. Some implications of these instantons to open inflation models are discussed in a separate paper 






























In addition, we have = 0 as mentioned at the end of the previous section. Hence, for arbitrary 64> of 0(e) we 

have 


+ 5(j)] 


5<t) 




Then the derivatives of S[(j)(; + Scj)] with respect to ( are expressed as 


dC 


d(j)Q 


d^S d(j)(^ d4)(^ 

Ic ^ ~W~ 


6(i) 




’ dC" . ■ 


(4.7) 


(4.8) 

(4.9) 


Now we analyze the behavior of the action under a variation of C and its relation to the value of the lowest eigenvalue 
for the C)(4)-symmetric fluctuations. 

First we show that dS[(j>i^ + 54)]/dC, = 0 if + d;/) is a solution to the equation of motion to 0{e). Note that a 
bounce solution generally cease to exist for continuous values of but it exists only at discrete values of C. Assuming 
that (j)c^ + (50 is a solution, we have 




Then dS/dC, in Eq. ([4.8|) can be expressed by using Eq. (4.10) as 


(4.10) 




525 ( 0 ) 

(502 


\<t>0 


(4.11) 


and this vanishes because of Eq. (4.5). 

Next we consider the shift of the lowest eigenvalue, A, due to and show that d?S/dC,'^ has the same signature 
as A does. In the unperturbed case, the lowest eigenvalue is zero, A = A^°^ = 0, whose eigenfunction is given by 
(^( 0 ) = d 0 ^/d^. To evaluate the shift of A, we consider the eigenvalue problem given by 


5^S 

6(j}^ 


[0c + <50;*,<7’] 


Xal (fi. 


(4.12) 


In particular, we have 


^[0f + (50; ip,ip\=xj al \ip\^dr]. 


(4.13) 


We denote the 0(e) correction to ip by (p = Using Eq. (^), the left hand side of Eq. ( 4.13 ) reduces 

to 




(4.14) 


The first term in the right hand side of the above equation can be rewritten by using Eqs. (4.6) and (4.10) as 
535 ( 0 ) 


(503 

Thus, to 0(e), we obtain 




r 52 q(i) 

A J al\ip\^df] =+ 


525(0) 

d02 


0^;50, 




dC 


6(j) 


4’(i 


Scj) 

(i(p(3) 


0 c; 




dC 


(4.15) 


dC 


d^S 

dC 


(4.16) 


where the last equality follows from Eq. (4.9). Therefore the sign of the eigenvalue A is the same as that of d^S/dC,^. 
This tells us that when there appears a negative mode, i.e., if A < 0, the action of the bounce solution is not a local 
minimum. 
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To gain a more definite picture of the behavior of the action, we consider possible models for with a certain 
variable parameter. Among all the possibilities, we give schematic plots of three typical cases of + S(j)] = [(/)<;] 

in Fig. ||(a) - (c). In each plot, the dotted line represents the case of a parameter for which the bounce solution has a 
positive (PS/d(^^. As the parameter is varied, the shape of smoothly passes through the one like the dashed 

line and changes to the one given by the rigid line. In the case (a), the value of ^ at which takes an extremum 

value shifts toward zero, and only the Hawking-Moss type bounce at C = 0 is left in the end. In the case (b), the 
value of C at the extremum of 5'^^^ [</)(;] does not shift much, but the amplitude of changes. At the moment 

when 5’^^^ is given by the dashed line, infinite number of solutions become degenerate, just like the case of the 
unperturbed action Beyond this critical point, there appears again one non-trivial solution, but with negative 

(PS/dC,^. The action is manifestly greater than that of the Hawking-Moss solution at C = 0. In the case (c), at the 
moment when is given by the dashed line, there appears a A = 0 mode at the extremum of 5"^^^ Beyond 

this critical point, there appear three non-trivial extrema. Two of them that newly appeared have positive values 
of d^S/dC,"^ while the one corresponding to the original solution has negative d‘^S/d(‘^, hence has a negative mode. 
Also in this case, the latter solution with a negative mode has a larger action than that for either of the other two 
solutions. 



FIG. 3. Three typical cases of the behavior of the Euclidean action as a function of when a model parameter is varied. 


As an example of 5”^^^ let us consider the one given by 


5(1) = e 


4 

dry Vi, Vi= a„ 

n=0 


(4.17) 


where 


2{(i)/4>^) - 1 . 

5(i)[</)c] = 


In this case, A^i) [0^] can be calculated analytically as 
e |(ao + a4) + ^^^T-T^ 3 ^(coshCsinhC-C) + ^^ 

.3 smh C smhC 


(( coth ( 



(4.18) 


If we choose a one-parameter family of models, say, like 2 (q;i -I- aa) = —2 -|- a and 4a2 = 3, where a is the variable 
parameter, the situation as shown in Fi^H(a) is realized. On the other hand, if we choose 2(ai -I- era) = —a and 
4a2 = 2a, the situation as shown in Fig. Mb) is realized. Within this limited class of perturbations, the situation as 
shown in Fig. |^(c) is not realized. 


V. CONCLUSION 

We have investigated the negative mode problem associated with false vacuum decay with gravity. We have shown 
that the existence of a negative mode around a non-trivial Euclidean solution, called the Coleman-De Luccia bounce 
solution, is equivalent to that of a supercritical supercurvature mode in the perturbation spectrum of the quantum 
fluctuations in the open universe that appears inside the bubble. Supercritical supercurvature modes are those for 
which the mode functions diverge exponentially for large spatial radius in the open universe. Then we have proposed 
a conjecture that there exists no supercritical supercurvature mode. If this is true, there will be no negative mode 
around the Coleman-De Luccia bounce solution that dominates the process of false vacuum decay. 

To investigate the validity of our conjecture, we have first provided a potential model that admits a continuous 
series of Coleman-De Luccia type bounce solutions under the weak gravitational backreaction approximation. This 
series contains a Hawking-Moss type solution as a limiting case and each of these bounce solutions has a zero mode 
as the lowest eigenvalue. Then for a class of potentials that can be realized by small modifications of this potential 
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model, we have analyzed the behavior of the Euclidean action around a bounce solution. We have shown that its 
Euclidean action is not the smallest among non-trivial solutions if there exists a negative mode. 

We have considered three typical cases of the behavior of the action when a model parameter is varied. In all of 
these cases, we have found that, when there appears a bounce solution with a negative mode, it does not give the 
smallest action but there exists another bounce solution, either Coleman-De Luccia type or Hawking-Moss type, with 
a lower action that has no negative mode. This evidence strongly supports the conjecture that there is no negative 
mode for the Coleman-De Luccia bounce solution that dominates the tunneling process. 
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